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SUMMARY / " / 

This DaDer explains how NASTRAN can be utilized advantageously in the preliminary design cycle. The 
ink™ Sn of Z preliminary design process lends itself to programs that can produce multiple 
configurations or variations on a particular design with minimal cost or effort. The latter portion of the 
nr^esHncSmoassTs refining the design and adding more detailed analyses (particularly for other 
disciplLs). A method for quickly generating balanced spacecraft loading conditions for use in preliminary 

design and analysis also is explained. 

The following additional sections are included: 


1. Background 

2. Symbols 

3. Analytical Process 

4 . Aerodynamic Load Distributions 

5. NASTRAN Applications 

6. Conclusion 

7. References 


BACKGROUND 

The preliminary design cycle seeks to obtain general as well as specific information rapidly and 
inexpensively ye/accurately. The preliminary design cycle (see fig. 1) for spacecraft or space systems 

usually involves evaluating multiple designs for a given configuration or ^Zused 

configurations. A process for the analysis and evaluation work has been established (ret. l) ana usea 
(ref 2) for several investigations. This process (fig. 1) starts with a solid representation of the design and 

evoives into a finite element representation for static and dynamic ^^^^VtraN The p^essrf 
nerforminc the finite element analysis. Two such systems are IDEAS and NASTRAN the process oi 
preliminary design has, among other things, two objectives that can be opposing: (1) to provide an 
analytical representation that can be easily revised, and (2) to provide ™ 

refined as pan of the design improvement after a configuration has been accepted. The IDEAS system 
readHy knthf itself to objective number (1), while NASTRAN is part.cularly useful for objeettve 

number (2). 

Various researchers have suggested approaches (ref. 3 and 4) for optimizing a s^ctu^ design The 
optimization researchers usually stan with a given configuration and loading condition. The preliminary 
design issues addressed in this paper allow consideration for a broader viewpoint. This broader viewpoint 

asks the following questions: 


1 . What is a good configuration? 

2 . What vehicle loads go with a particular configuration? 
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SYMBOLS 


The following symbols are used in this paper: 

air density (slugs/ft^) 

angle of attack (rad) 
computer-aided design 
dynamic pressure (lb/ft^) 
gimbal angle (deg) 
gust velocity (ft/sec) 
normal force coefficient slope (1/deg) 

payload 

solid rocket motor 
surface reference area (ft^) 
thrust (lb) 

vehicle velocity (ft/sec) 
wind velocity (ft/sec) 

X coordinate of the center of gravity (in.) 

X coordinate of the center of pressure (in.) 

X coordinate of thrust vector application point (in.) 

ANALYTICAL PROCESS 

The process starts with a candidate design or configuration that needs to be evaluated. A computer-aided 
design (CAD) representation is created and serves as a basis for the finite element model. The basic finite 
element model can serve as the starting point for investigating alternate configurations. It usually takes at 
least one iteration through a segment (see fig. 2) of the process to get a reasonable estimate of the structural 
sizing and weights. The first pass-through also provides a good test of the model fineness. The analyst 
would like the finite element model to be fine enough to give believable stress and deflection predictions; 
however, it should be crude enough to keep computing costs and time at a low level. 

The box entitled “Finite Element Model” (see fig. 2) could utilize any one of a number of different 
programs. The two most attractive systems for this project were IDEAS and NASTRAN. Table I gives a 
comparison between the two systems. In order to generate a good preliminary design, both programs (or 
other comparable ones) should be used: IDEAS (to compare configurations and to select one) and 
NASTRAN (to provide the starting point for detailed design and certain specialized analyses [e.g., flight 
control, flexible body loads, etc.]). This is shown in fig. 3. 

Considerable effort has been spent in computing vehicle load conditions that are configuration 
dependent. Any preliminary design can only be as good as the vehicle loads being used. The issue of 
balanced load conditions is important because in the early stages (preliminary design) of a design 
meaningful loads are very difficult to obtain. Balanced load conditions on a vehicle allow an analyst to look 
at the computed stresses and deflections and not be concerned about how the results have been skewed by 
assumed boundary conditions or unbalanced loads. A balanced load condition is one where the sum of all 
forces and moments (aerodynamic, inertial, and thrust) acting on the vehicle are zero. 
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Table I. Comparison of IDEAS and NASTRAN 



FINITE ELEMENT MODEL 


IDEAS 

Advantages 


Quick turnaround 
Complete solutions 
Low cost 

Static and dynamic results 
Color graphics 
Database capability 

iOtsadvantages 

Limited capability to interface with other 

programs/disciplines 

Limited usage in U.S. 

Available on limited platforms 

Specialty (e.g., buckling, etc.) solutions not available 


Easy to interface with other programs/disciplines 
Wide usage in U.S. 

Highly portable 

Sophisticated solutions available 


Not easy to generate multiple configurations 

No built-in color graphics 

No convenient database features 

'Not so quick’ turnaround 

Not particularly low cost 


AERODYNAMIC LOAD DISTRIBUTIONS 

An auxiliary program was set up to P™ vide A'** 1 ' for evaluating 

These uncertainties may ^^ctory 

parameters, such as: dynamic pressure (q), angle of attack (a), or vehicle center of gravtty locatton <Xcg>- 

Aerodynamic forces normal to the vehicle lonp. udina, ■ ^ ^""o "Te 

the vehicle structure. They also requtre <^™^'/Xle As^hown in fig 4 the loads analysis uses 

Aerodynamic inputs consist of the normal md^mStrond on the size'amd sha^of *e vehicle 

pressure (q). The vehicle size and shape % Q g is Really represented by distributed normal 

CN a distributions are oh^ned empirically or from tesjdau 

available for similar configurations Empirical methods (ref. 5) were used for *' C N« 
variations along vehicle components of various shapes and for a wide range of flight Mach numbe . 

The magnitude of a is typically obtained from dynamic trajectory simulations with - 
shear and gusts. If trajectory simulations are no. " “ ^pproxtmate value for a can be 
superimposing the wind and gust speeds (ref. 6) on the vehicle speed. 

a = T an-l(- ^ n ^ V g^ \rad (1 > 

V V vehicle ) 
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With the trajectory parameters of q and a, and with the Cjq a distribution defined along the vehicle, the 

auxiliary program is used (fig. 5). The method computes the distributed normal forces, net pressures, and 
the summed forces and moments about the vehicle's center of gravity. Using this method, a vehicle segment 
of incremental length is subjected to an aerodynamic normal force where the magnitude depends on CN a , 

q, and a. 

ANormal Force = q S re f Cjsj a a, lb/in. (2) 


where: 


Cjq a = distributed normal force coefficient slope, l/(in.-rad) 

q = dynamic pressure, 1/2 p (V vc hicle) 2 Ob/ft 2 ) 

S re f = reference area (ft 2 ) 

a = angle of attack (rad) 
p = atmospheric density (slugs/ft 2 ) 

The above equations are used to compute the normal load distribution along the vehicle. It is then integrated 
within the auxiliary program to compute the load and moment summations about the center of gravity. The 
presence of additional elements, such as solid rocket motors (SRMs), can be accounted for by adding their 
point-load contributions to the total forces and moments. 

Static balance calculations are included in the program to determine the amount of engine gimbal angle 

(5) required to overcome (or balance) the aerodynamic moment. This is computed from the moment balance 
between the aerodynamic forces and the engine thrust, as shown below. 

TSin(S)(Xgj m t> a i “ X C g ) = X(CN a )q«Sref (^cg ~ ^cp ) (3) 


The above equation is then solved for the gimbal angle, 5. 


5 = Sin * 


(X C g — X cp ) 

T(Xgimbal _ ^cg) 


I- deg 


(4) 


where: 

I(C Na ) = integrated normal force coefficient slope on vehicle (rad) 


q = dynamic pressure (lb/ft 2 ) 
S re f = reference area (ft 2 ) 

T = engine thrust (lb) 

X C g = center of gravity station (in.) 
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X cp = center of pressure station (in.) 
Xgimbal = engine gimbal station (in.) 


a = angle of attack (rad) 

6 = engine gimbal angle for balancing the aero forces (deg) 

For the case when additional engines exist, as in the case of SRMs, the above static moment balance is 
altered to include such engines. With the SRM and Core subscripts used for the appropriate elements, the 
moment balance expression becomes: 


OCore + t SRM )S in(6)(X g j m bal - X cg ) - qaS rc f { c N aCore ( x cg “ X cpCore 

+Cn «SRM (X cg _Xc PSRM } 


( 5 ) 


where Cn corresponds to the core stage element and is equivalent to Z(Cw ) in the previous 

a Core 

moment balance equation. 

Then: 


5 = Sin * 


qaS r ef{cjM 


aCore (Xc S ~ Xc PCore } + CN asRM (Xcs ~ Xg PSRM } 
^ Core ( x gimbal ~ x cg)'*'^'SRM( x gimbal — x cg) 




( 6 ) 


With the gimbal angle defined, the axial and tangential thrust values are calculated. These thrust 
components are then used to compute the axial and tangential accelerations (normal to the vehicle 
longitudinal axis), which are input into the finite element model. 

, total axial thrust 

axml acceleration = vehk]e welght 

, total tangential thrust + Z( normal force) /os 

tangential acceleration = vehicle weight <8> 

Key load parameters can be changed easily in the program to see their influence on loads and engine 

control deflections. A change in dynamic pressure, (q), angle of attack (a), or vehicle center of gravity 
(Xc g ) will readily show the sensitivity of aerodynamic loads to such changes. 

NASTRAN APPLICATIONS 


The Background section of this paper discussed using two different finite element programs for 
structural analysis. Why not just use one model/program for the entire preliminary design cycle? The two 
systems, IDEAS and NASTRAN, have different advantages and disadvantages (see table I). 

The finite element solver that is internal to IDEAS is a valuable tool, especially when rapid results based 
on model variations are desired; however, for certain applications, a NASTRAN finite element 


103 



representation is much more useful. Figure 6 shows some of the static and dynamic applications that can be 
supported by the NASTRAN model. 

The IDEAS finite element model can be used in its full mass and stiffness representation to compute the 
first few system mode shapes and natural frequencies of the accepted configuration. This information can be 
used as a check on the mode shapes and frequencies that are later computed using a reduced dynamic model 
(e.g., flexible body loads model) generated with NASTRAN. 

CONCLUSION 

The early portion of the preliminary design cycle makes the use of the finite element code in IDEAS 
attractive because a vehicle analysis can be quickly redone after sizing changes are made. This paper 
describes a procedure for preliminary design and shows how NASTRAN can be used as a vital tool in that 
process. Additionally, a method for setting up balanced vehicle load conditions, as an integral part of that 
procedure, has been explained in detail. The challenge in the preliminary design cycle is to create a large 
amount of meaningful information rapidly and inexpensively, to use the preliminary design analytical 
representation to interact with many disciplines, and to support the evolution of a detailed design. 

The later stages of the preliminary design can be effectively handled by NASTRAN because of its ability 
to: 

1 . Handle many thousands of degree problems relatively cheaply 

2 . Run on many different platforms 

3 . Easily interface with other programs/data sources 
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Figure!. Preliminary Design and Analysis Cycle 
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Figure 2. Units Loads are Used to Assess the Initial Design 
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Figure 3. NASTRAN and IDEAS are Used for Different Purposes 
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Figure 4. An Example of High Dynamic Pressure Region Vehicle Loads 
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Figure 5. The Aerodynamic Influence is Displayed Different Ways 
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